Peculiar scaling of self-avoiding walk contacts 
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The nearest neighbor contacts between the two halves of an iV-site lattice self-avoiding walk offer 
an unusual example of scaling random geometry: for N — > oo they are strictly finite in number 
but their radius of gyration R c is power law distributed cx R^ T , where r > 1 is a novel exponent 
characterizing universal behavior. A continuum of diverging lengths scales is associated to the R c 
distribution. A possibly super-universal r = 2 is also expected for the contacts of a self-avoiding or 
random walk with a confining wall. 
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The self-avoiding walk (SAW) is a classical problem in 
statistical mechanics, playing a central role in our under- 
standing of polymer statistics and intimately related to 
magnetic critical phenomena and percolation jjj. In its 
most simple version the SAW amounts to the statistical 
characterization of equally probable single chain confor- 
mations, with no overlaps or intersections. These con- 
formations are made by N — 1 successive nearest neigh- 
bor (nn) steps (N sites) on a lattice in d dimensions. 
For N — > oo, quantities like the radius of gyration with 
respect to the center of mass of SAW configurations, 
R g , have an average (R g ) ~ i\T I/SAW , where ^ SAW is the 
SAW metric exponent. When d > 1 self- avoidance does 
not prevent the conformations from involving close ap- 
proaches of different parts of the chain: here we generally 
count as contacts the pairs of nn lattice sites which are 
visited, not consecutively, by the SAW. The totality of 
such contacts grows on average proportional to N and 
possesses the same fractal dimension (= l/u SAW ) as the 
whole SAW. 

In the present Letter we discuss so far unexplored fea- 
tures of a particular subset of SAW contacts. Such fea- 
tures represent an unusual example of how scale invari- 
ance can manifest itself in a finite, non-extensive por- 
tion of an infinite fractal set. The scaling exponent r 
of the probability distribution of the gyration radius of 
such subset represents a novel characterization of SAW 
universal behavior. The fact that r > 1 implies the ex- 
istence of a whole continuum of diverging characteristic 
lengths in the SAW problem, in addition to the length 

(Rg)- 

The subset on which we focus here is that of the con- 
tacts between the two halves of a SAW (Fig. [j]). Count- 
ing such contacts alone allows to get rid of less interesting 
effects, which are extensive in N . Being related to prob- 
lems like network formation, transport or intramolecular 
reactions, these contacts can be of particular interest for 
applications in which the two half-chains are made of dif- 
ferent monomers, like for diblock copolymers [Q. When 
the chain represents a homopolymer, the subset we con- 
sider is particularly significant in relation to the effect 



of nearest neighbor interactions on the SAW jj). This 
is the case of models of the polymer collapse |l],|||3), 
where an attractive nn energy e < is associated to each 
contact. In this case a Boltzmann factor e _£ / T weighs 
each contact occurring in a configuration at temperature 
T. If for such model one counts only contacts between 
the two halves of the SAW, the average number of them 
is of the order in the high T regime, increases as 
with < (j) < 1, at the point, and scales as N at low 
T. <t> turns out to coincide with the crossover exponent 
e at T = T e , where (R g ) ~ N Ue , with v e £ v SKW [§. 
Similar behaviors occur if the attractive interactions are 
associated exclusively with this subset of contacts, and 
the model describes a diblock copolymer zipping transi- 
tion @,|. In different T regimes of the O and zipping 
transitions the contacts between the two halves of the 
SAW behave in a way analogous to the monomers adher- 
ing to a wall in polymer adsorption J7| . Most recently it 
has also been found that in d = 2 the fractal dimensions 
of the contacts between the two half-chains are the same 
at the homopolymer point and at the diblock copoly- 
mer zipping transition Q. Thus, intriguing universality 
aspects can be expected to underlie the geometry of SAW 
contacts. 

The contacts between the two halves of a polymer have 
already been studied by renormalization group meth- 
ods H in the high T regime controlled by excluded vol- 
ume. The focus there was on the scaling of their aver- 
age number, (N c ), and precisely on the scaling correc- 
tion exponent describing its approach to a finite limit for 
N — > oo. This finite limit shows that only a vanishing 
fraction of the total number of contacts pertains to ap- 
proaches of remote portions of the chain. In fact (N c ) 
gives only limited information on the contact statistics. 
The full probability distribution function (PDF) for N c , 
P(N C , N), could in principle have higher moments di- 
verging with TV — > oo, and we devote a first effort to 
check this possibility, which is normally not considered 
in polymer statistics. We perform several Monte Carlo 
simulations to investigate this PDF and its moments in 
various dimensions and for 80 < N < 2000. Since the 



1 



configurations with many contacts are very rare, we em- 
ploy a pruned-enriched Rosembluth method (PERM ||) 
tuned to increase the sampling of configurations with 
high number of contacts. This enables us to obtain a 
good statistics also in the tail of the PDF, which for all d 
turns out to be a negative exponential without substan- 
tial dependence on N, as displayed in Fig. || for SAW on 
cubic lattice. 

These results show that this set of contacts is strictly 
finite in the N — > oo limit. In spite of this, one can 
still ask whether these contacts have interesting geo- 
metrical scaling properties, not just reducing to those 
of a spatially bounded random set. If we indicate by 
Prad(Rc,N) the cumulative PDF of R c (Fig. |l|) over all 
SAW configurations, strict boundedness would mean that 
the moments (Rf) ~ J dR c R1P ra d(R c , N) remain finite, 
for N — > oo, Vg. We extrapolate the moments in the 
form (R q ) ~ N" q . The data are generated by sampling 
SAW of fixed length with a Monte Carlo algorithm based 
on pivot moves [p.0| , which have been proved to be very 
efficient for SAWs Since the contacts between the 
two halves of the SAW are mainly located close to the 
junction point, local moves [ fl2"| are also attempted in this 
region. We consider hypercubic lattices and the FCC lat- 
tice in d = 3. Contrary to what happens for P(N C ,N), 
for all d we find that a q is positive and grows linearly with 
q, at sufficiently high q. The data are consistent with a 
behavior 

JO if q<r-l m 

Gq ~ \ v\q-{T-\)\ if q>r-l , W 

where a q — may represent logarithmic divergences for 
< q < t — 1. In d = 2, for example (see Fig. ||), 
we find t = 1.93(2) and v — 0.76(1), consistent with 
3/4 jn|. The form @) is compatible with a 
PDF having a scaling form P ra d{Rc, N) ~ R~ T f{R c /N u ). 
That v — i/ SAW is quite plausible since jV" SAW is the 
only expected characteristic length in the problem. How- 
ever, as discussed below, P ra d itself introduces a mul- 
tiplicity of new lengths scales for the SAW. The fact 
that the moments do not approach zero for q < r — 1 
(i.e. u q does not become negative) should be due to the 
circumstance that /(x) is not converging to zero for 
x — > 0. In other terms, the scaling function g, if we 
write P ra d — N~ T "g(R c /N l/ ), is singular for its argu- 
ment approaching zero (g(x) ~ x~ r ). We verify these 
assumptions on the structure of P ra d by scaling collapse 
plots for various d. For example, Fig. |j shows the data 
collapse of ln[/(x)] for d — 3. A similar collapse plot for 
ln[g(x)] is reported in the inset. 

It is indeed the singular character of g(x) for x — > 0, 
which allows the exponent r to take a nontrivial value 
> 1, while maintaining the zero-th, normalization mo- 
ment of the PDF equal to 1. The lower length cutoff I 
(lattice spacing in this case) is crucial to obtain a finite 
integral, in dx = d(R c /N u ), of the PDF in the contin- 
uum limit, because the main contribution comes from 



small values of x, close to x_ = l/N w . This contribution 
has an iV-dependence which compensates the diverging 
factor N~ V ( T ~ X } extracted in front of the integral. 

r > 1 means that we can not associate to the R c 
PDF a unique characteristic length. Indeed, putting 
= (R q c ) 1/q we find £ 9 - N v * with v q = a q /q, for 
q e [t — l,+oo). This means that the self-similarity of 
contacts has an intrinsic multiscaling character. 

The r found here is a novel exponent for the SAW. It 
is a measure of the spread of the region within which one 
half of the chain feels the presence of the other one in 
the SAW configurations. A higher r (see Table |) indi- 
cates more localized contacts. Like the global geometry 
of the SAW, the spread of contacts is determined by the 
interplay between entropic and excluded volume effects. 
It is remarkable the non-monotonic behavior of r, which 
takes minimum values in d = 3 and d = 4, indicating 
these dimensionalities as the optimal ones for a broad 
interpenetration of the two half-chains. At high d there 
is soon much room for the two branches of the SAW to 
develop without giving rise to contacts, and this tends 
to localize them more, r is rather high in d = 2: we in- 
terpret this as a consequence of the peculiar topology of 
the two-dimensional lattice, which makes it more difficult 
than, e.g., in d = 3 for the two half-chains to approach 
each other at large length scales. The dependence of r 
on d, also above the upper critical dimension d u = 4, is a 
further indication of the peculiar novel character of this 
exponent. 

In d=2 and d=3 we also investigate the behavior of 
the contacts between the two half-chains in the presence 
of nn attractive interactions (8 point model). While for 
T > T e the behaviors of their PDF's appear the same as 
at T = oo, at the point the moments of P diverge as 
those of P ra d, while r becomes equal to 1 and <r q / q = v e 
for the latter PDF. This suggests that the disappear- 
ance of the singular scaling function in P ra d could be a 
good criterion for locating the transition point. This is 
illustrated in Fig. [|, referring to the O point in d = 3, 
simulated as in ref. || with chains up to N — 2000. We 
collect data from two runs at —e/T = 0.25, 0.27 and we 
use the multiple histogram method [|l4| to calculate the 
moments in a surrounding interval of temperatures. The 
result from the values of q examined is —e/T = 0.274(4), 
consistent with accurate estimates by other methods: for 
example, 0.275(8) in § and 0.2690(3) in §. Similar 
results are valid for models of the diblock copolymer zip- 
ping transition |jq| . 

For SAW with attractive interactions the average num- 
ber of contacts between the two half chains behaves in the 
various temperature ranges as the mean number of SAW- 
wall contacts in a polymer adsorption transition. Thus, 
it makes sense to check whether also in the adsorption 
process the high T ordinary regime is characterized by 
peculiar scaling features of the kind described above. To 
this purpose we study the PDF's of the number and the 
radius of SAW- wall contacts at T — oo, where the wall 
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exerts only a geometrical confinement effect on the chain. 
Here we consider as radius the mean distance R* of a 
contact from the point where the SAW is grafted to the 
wall J16j . Also in this case all the moments of the PDF 
of N c appear to remain hnite for N — > oo. The PDF's 
of R c show singular scaling functions with an exponent 
t ~ 2, in d = 2 and d = 3 (Table ||). 

We study also the case of RW confined by a wall, for 
which we are able to compute exactly r. Consider first 
a RW on a square lattice tilted of 45° with respect to 
the coordinate axes, in such a way that the nn of a site 
have coordinates {(1,1), (—1,1), (—1,-1), (1,-1)}. In 
this way a nn step moves the RW with nonzero and inde- 
pendent displacements along both coordinate directions. 
This simplifies the calculations, but we expect the fi- 
nal results to be valid for every lattice model, because 
independence is asymptotically recovered for long RW. 
The generalization to d dimensions gives 2 d nn vectors of 
the form (1, 1, . . . , 1), (-1, 1, . . . , 1), . . . , (-1,-1,..., -1) 
and, for example, in d = 3 one obtains the BCC lattice. 
Let the walk start from the origin, near a d — 1 dimen- 
sional hard wall perpendicular to the x coordinate: the 
problem in the x direction is equivalent to a one dimen- 
sional RW which steps to nn sites with equal probability. 
The probability to be again at x = after n steps is 
given by P„(0) = 2""( " ) . The effect of the impen- 
etrable wall is represented by forbidden sites located at 
x = — 1. So one has to subtract the probability to travel 
through this point. The method of images |17| ] in this 
simple case says that this is equal to the probability to 
go from x — to its mirror image with respect to x, 
x = -2. So P„(0) = P„(0) - P„(-2) = P„(0)2(n + 2)" 1 
is the probability to be on the surface. The q-th. moment 
of the average distance of a contact from the origin is 
thus 
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where n 1 / 2 is the root mean square displacement after 
2tt n n+1 / 2 e~ n for the binomials, 



n steps. Using n! f 
one recovers that the denominator is finite for N — > oo, 
while the numerator is equivalent to a sum of the type 
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diverging as A^s- 1 )/ 2 if g > 1. This 
means r = 2 for RW near a wall, in any d. The ex- 
act result for RW suggests that, in any dimension, for 
SAW confined by a wall, the scaling of the contacts is 
the same as for RW. Thus, excluded volume effects seem 
to play no role in determining r for the contacts of the 
SAW with a confining wall. 

In summary, the contacts discussed above represent a 
very peculiar example of scaling random set: in fractal 
physics we are familiar with sets in which the number of 
elements is growing to infinity together with their aver- 
age radius of gyration, and criticality implies a nontrivial 
scaling for the PDF's of both quantities. A typical ex- 
ample are percolation clusters Q , whose size PDF has a 



singular scaling function with nontrivial r, like we instead 
find here for P ra d- In the case considered here the scale 
invariance of the set is not accompanied by its number 
of elements being broadly, power law distributed. The 
criticality is in fact triggered by the length of the whole 
chain, N, becoming infinite, while N c remains finite. A 
t exponent, possibly super-universal, can also be defined 
for the contacts between a SAW and a d — 1 dimensional 
confining wall in the ordinary regime. The case of a RW 
in presence of wall can also be treated, yielding the exact 
classical value r = 2, which could apply also to SAW, 
independent of d. Exact determinations of these expo- 
nents for SAW in low (d = 2) or high dimensionality, and 
the possible descriptions of the new scalings within the 
renormalization group framework remain a challenge for 
the future. 
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TABLE I. Extrapolated values of r and v for various lat- 
tices, using SAW with length from N M in to N M ax- 
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TABLE II. As in table @, but for SAW- wall contacts 




FIG. 1. SAW in two dimensions: the contacts between its 
two halves (light and dark, respectively) are indicated by open 
circles. CM. is the center of mass of these contacts and R c 
their radius of gyration. 
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FIG. 3. (j q vs q and extrapolation of r — 1 for d — 2. We 
expect o q — for q < t — 1. Numerically a logarithmic diver- 
gence can not be easily distinguished from a power law one 
with a q > 0. 
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FIG. 4. Collapses of ln[/(a;)] (see the text) and ln[gr(a;)] (in- 
set) for SAW on cubic lattice, with N = 200 (O), 400 (X), 
800 (A) and 1600 (□). 



0.65 




t;so°°°" | 

FIG. 2. Histograms of \n[P{N c , N)] for SAW on a cubic 0.24 0.25 0.26 0.27 0.28 0.29 

lattice, with N = 200 (O), 400 (X), 800 (A) and 1600 (□). ~ e/T 

FIG. 5. Extrapolation of <Jq/q close to Te in d — 3, for 
q = 0.5 (O), Q = 1 (X), q = 2 (A) and q = 4 (□). The 
crossings are consistent with the expectation to find r = 1 
and Oq/q = ^ e = 1/2 right at T — T e - The curve with the 
opposite trend (+) refers to the effective v exponent of the 
SAW radius of gyration. 
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